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Abstract. Calculations of nonequilibrium processes become increasingly 
feasable in quantum field theory from first principles. There has been impor¬ 
tant progress in our analytical understanding based on 2PI generating function¬ 
als. In addition, for the first time direct lattice simulations based on stochas¬ 
tic quantization techniques have been achieved. The quantitative descriptions 
of characteristic far-from-equilibrium time scales and thermal equilibration in 
quantum field theory point out new phenomena such as prethermalization. 
They determine the range of validity of standard transport or semi-classical 
approaches, on which most of our ideas about nonequilibrium dynamics were 
based so far. These are crucial ingredients to understand important phenomena 
in high-energy physics related to collision experiments of heavy nuclei, early 
universe cosmology and complex many-body systems. 
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1. Characteristic far-from-equilibrium time scales 

Understanding the dynamics of quantum fields far away from the ground state or 
thermal equilibrium is a challenge touching many aspects of physics, ranging from 
early cosmology or collision experiments with heavy nuclei to ultracold quantum 
gases in the laboratory. One of the most crucial aspects concerns the characteristic 
time scales on which thermal equilibrium is approached. Much of the recent interest 
derives from observations in collision experiments of heavy nuclei at RHIC. The 
experiments seem to indicate the early validity of hydrodynamics after somewhat 
less than 1 fm, whereas the present theoretical understanding of QCD suggests a 
longer thermal equilibration time. 

Here it is important to note that different quantities effectively thermalize on 
different time scales and a complete thermalization of all quantities may not be 
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Fig. 1. Fermion occupation number for three different momentum modes as a 
function of time in the chiral quark meson model of Ref. [ CP ■ 

necessary to explain the observations. This has been pointed out in Ref. [[J, 
where it was shown for a chiral quark-meson model that the prethermalization of 
important observables occurs on time scales dramatically shorter than the thermal 
equilibration time. As an example, Fig. 0 shows the nonequilibrium time evolution 
of fermion occupation number for three different momentum modes in this model. 
The evolution is given for two different initial particle number distributions A and 
B shown in the insets, with same energy density. Therefore, both runs have to 
lead to the same distributions in thermal equilibrium. The vertical line in Fig. ^ 
marks the characteristic time scale ~ idamp, after which the details about the initial 
distributions A or B are effectively lost. One observes that this happens far before 
the late-time approach to thermal equilibrium. The time tdamp for the effective loss 
of initial condition details is an important characteristic time scale in nonequilibrium 
dynamics, which is very different from the thermal equilibration time t eq . Note that 
the long-time behavior to thermal equilibrium is shown on a logarithmic scale in 
Fig. H ( in units of the scalar thermal mass to). 

In contrast to the very long time t eq for complete thermal equilibration, prether¬ 
malization of the (average) equation of state sets in extremely rapidly on a time 
scale 

^pt ^ fdamp ^ teq ■ (1) 

In Fig. Owe show the ratio of average pressure (trace over space-like components of 
the energy-momentum tensor) over energy density, w = p/e, as a function of time. 
One observes that an almost time-independent equation of state builds up very early, 
even though the system is still far from equilibrium! Here the prethermalization 
time t p t is of the order of the characteristic inverse mass scale to - 1 . This is a typical 
consequence of the loss of phase information by summing over oscillating functions 
with a sufficiently dense frequency spectrum. If the “temperature” (T), i.e. average 
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Fig. 2. The ratio of average pressure over energy density w as a function of time. 
The inset shows the early stages for two different couplings h and demonstrates that 
the prethermalization time is rather independent of the interaction details. 

kinetic energy per mode, sets the relevant scale one finds T i pt ~ 2 — 2.5 [HJ. For 
T > 400 — 500 MeV one obtains a very short prethermalization time t pt of somewhat 
less than 1 fm. 

This is consistent with very early hydrodynamic behavior, however, it is not 
sufficient as noted in Refs. [□E] Beyond the average equation of state, a crucial 
ingredient for the applicability of hydrodynamics for collision experiments [Ej is the 
approximate isotropy of the local pressure. More precisely, the diagonal (space-like) 
components of the local energy-momentum tensor have to be approximately equal. 
Of particular importance is the possible isotropization far from equilibrium. The 
relevant time scale for the early validity of hydrodynamics could then be set by the 
isotropization time. The analysis of simple models lead to an isotropization time 
which is given by the characteristic damping time ~ tdamp (cf. Figs. El [EJ- However, 
a possible weak-coupling mechanism for faster isotropization in gauge theories such 
as QCD has been identified in Ref. [0 in terms of plasma instabilities [EJ. Whether 
this can explain the experimental observations or whether they suggest that we 
have to deal with some new form of a “strongly coupled Quark Gluon Plasma” is 
an important open question. 

2. What can we learn from transport or kinetic theory? 

For out-of-equilibrium calculations there are additional complications, which do not 
appear in thermal equilibrium or vacuum. Standard approximation techniques, such 
as perturbation theory, are not uniform in time and fail to describe thermalization. 
Aspects of systems with high occupation numbers are often successfully described 
using classical field theory methods. However, classical Rayleigh-Jeans divergences 
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Fig. 3. Snapshots at times t ~ 0, tdamp/2, fdamp and 3tdamp/2 of the initially 
anisotropic occupation number distribution as a function of p± (vertical) and p || 
(horizontal). The shown resolution was achieved by using a 64 3 spatial lattice. 

and the lack of genuine quantum effects — such as the approach to quantum thermal 
equilibrium characterized by Bose-Einstein or Fermi-Dirac statistics - limit their 
use. 

Most theoretical approaches to the important question of thermalization have 
been limited to semi-classical systems in the weak-coupling limit so far. With the 
advent of new computational techniques a more direct account of quantum field 
degrees of freedom becomes more and more possible. There has been important 
progress in our understanding of nonequilibrium quantum fields using suitable re¬ 
summation techniques based on 2PI generating functionals [ IJ. They have led to 
quantitative descriptions of far-from-equilibrium dynamics and thermalization in a 
variety of scalar and fermionic quantum field theories so far. An example for the 
results of such a first-principles calculation within a quark-meson model in 3+1 di¬ 
mensions is given in Figs. Q] and El One important application of these quantum field 
theoretical 2PI methods is to test standard transport or semi-classical approaches, 
on which most of our ideas about nonequilibrium dynamics are based so far. 

Several of these tests have been done for simple scalar A<f > 4 quantum field theo¬ 
ries for not too strong couplings A, which are well under quantitative control using 
2PI techniques. Following Ref. [0J, we consider a class of anisotropic initial con¬ 
ditions with an initially high occupation number of modes moving in a narrow 
momentum range around the “beam direction” p 3 = p y = ±p ts ■ The spatially 
homogeneous occupation numbers for modes with momenta perpendicular to this 
direction, p\ +p\ = p\, are small or vanishing. The situation is reminiscent of some 
aspects of the anisotropic initial stage in the central region of two colliding wave 
packets . 1 Of course, a peaked initial particle number distribution is not very specific 
and is thought to exhibit characteristic properties of nonequilibrium dynamics for 
a large variety of physical situations. 

An example of the earlier stages of such an evolution is shown in Fig.|3 Shown 
are snapshots of the occupation number distribution as a function of p± (vertical) 
andp|| (horizontal), where bright (dark) regions correspond to high (low) occupation 
numbers. The initial Gaussian distribution is centered around a momentum of the 

1 Ollier interesting scenarios include “color-glass”-type initial conditions with distributions ~ 
exp(— 1 p\/Qs) peaked around p 3 = 0 with “saturation” momentum Q s . 
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order of the renormalized thermal mass, and we will consider amplitudes of order 
one for the comparisons below (see [© for details). From Fig. [Hit can be seen that 
after the characteristic damping time fdamp the distribution starts to become rather 
independent of the momentum direction. Finally, at about 3f<jamp/2 the figure shows 
an almost perfectly isotropic situation. We emphasize that the distribution is still 
far from equilibrium. The situation is similar to what is displayed in Fig.^ where 
the vertical line indicates about fdamp- In general, we find that the characteristic 
damping time is well described by the standard relaxation-time approximation 

2d/ eq ) 47TTO ^ 

idamp - ~( eq ) — 0X2^2 ’ 

^Q 


where the imaginary part of the thermal equilibrium self-energy —E^ eq -*/2 in Fourier- 
space is evaluated for on-shell frequency o/ eq ) for momentum p ts . The second equal¬ 
ity in © is only valid for sufficiently high temperatures and weak couplings. We 
emphasize that the relaxation-time approximation © does not describe the ther- 
malization time t eq [2]. 

The question is whether transport or kinetic equations can be used to quan¬ 
titatively describe the early-time behavior t < tdamp, which is necessary if their 
application to the problem of fast apparent thermalization explained in Sec. Q is 
viable. The derivation of transport equations is standard. The evolution of particle 
number distributions is encoded in the time-ordered two-point correlation function. 
Its imaginary and real part is determined by the commutator and anti-commutator 
of two fields: 

p(x,y) =i([$(x),$(y)]), F(x, y) = i({$(z), $( 2 /)}) • (3) 

Here p{x, y ) denotes the spectral function and F{x, y) the statistical two-point func¬ 
tion. While the spectral function encodes the spectrum of the theory, the statistical 
propagator gives information about occupation numbers. Loosely speaking, the de¬ 
composition makes explicit what states are available and how they are occupied. 
For nonequilibrium F(x,y ) and p(x,y) are in general two independent two-point 
functions, whose exact time evolution equation reads 2 


[□* + M 2 (x)\ F(x,y) 



/ OO 

d 3 zT, p (x, z)F(z, y) 

-OO 

/ OO 

d 3 zY, F (x,z)p(z,y ), 

-OO 


(4) 


[□* + M 2 (x)\ p{x,y) = 



d 3 zT, p (x,z)p{z,y). 


(5) 


2 For a detailed derivation see e.g. Ref. [|H]. We are considering the equations for Gaussian 

initial conditions, which are underlying transport equations. More involved initial conditions can 
be considered using higher n-PI effective actions [151. 
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These are causal equations with characteristic “memory” integrals, which integrate 
over the time history of the evolution starting at time to = 0. Since they are ex¬ 
act they are equivalent to any kind of identity for the two-point functions such as 
Schwinger-Dyson/Kadanoff-Baym equations. Here Y,p{x,y) denotes the real part 
and —T, p (x,y)/2 the imaginary part of the self-energy E, where the local contri¬ 
bution is taken into account in M 2 {x). In 2PI approximations the self-energy E is 
obtained from the two-particle irreducible effective action [Gj. Here we consider a 
three-loop 2PI effective action m. It includes direct scatterings, off-shell and mem¬ 
ory effects. Most importantly in this context, it employs no derivative expansion. 
The latter is a basic ingredient for transport or kinetic theory. 

Transport equations are obtained from the exact equations by the following 
prescription: 

1) The lower bound (to = 0) of the time-integrals in ® is sent to the infinite 
past, i.e. to —> — oo. Of course, a system that thermalizes would have reached 
already equilibrium at any finite time if initialized in the remote past. Therefore, in 
practice a “hybrid” description is employed: transport equations are initialized by 
prescribing F , p and derivatives at a finite time using the equations with to —> — oo 
as an approximate description. 

2) Employ a gradient expansion. In practice, this expansion is carried out to 
lowest order (LO) or next-to-lowest order (NLO) in the number of derivatives with 
respect to the center coordinates X M = ( x p + y p )/2 and powers of the relative 
coordinates s p = x p — y p . 

3) Even for finite X° one assumes that the relative-time coordinate s° ranges 
from — oo to oo in order to achieve a convenient description in Wigner space, i.e. in 
Fourier space with respect to the relative coordinates. 

We emphasize that the ad hoc approximations 1) and 3) are in general not 
controlled by a small expansion parameter. They require a loss of information 
about the details of the initial state. More precisely, they can only be expected 
to be valid for sufficiently late times times t when initial-time correlations become 
negligible, i.e. (<f>(0, x)&(t, y)) ~ 0. The standard approximations 1) and 3) may, 
in principle, be evaded. However, if they are not applied then a gradient expansion 
would become too cumbersome to be of use in practical calculations. 

The NLO transport equations for F(X;p) = f d 4 s e lps F(X + s/2, X — s/2) and 
g(X;p) = —i J d 4 se ips p(X + s/2,X — s/2) in Wigner space are given by 

[2 p^d x , + (d x .M 2 (X))d Pll }F(X,p) = t e (X,p)F{X,p)-Z F (X,p)Q(X,p) 

+ { E F (X, p) , ReG fl (X, p) } pB + { ReE* (X, p ), F (X, p) } ^ , (6) 

[2p p d x » + (d x .M 2 (X)) d p „] e(X,p) 

= {£ e (X,p),ReG R (X,p)} pB + {Re£ R (X,p),g(X,p)} pB , (7) 

where G R (x,y) = 0(a; 0 — y°)p(x,y) and H R (x,y) = 0(a; 0 — y°)'E p (x,y) denote the 
retarded propagator and self-energy, and {.. ,} R b are the Poisson brackets. The 
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Fig. 4. Comparison of LO, NLO and full result for the on-shell 2 uid x oF with 
A = 0.5. The LO transport equation fails to describe the full results until rather late 
times. Taking into account the substantial NLO corrections the gradient expansion 
becomes quite accurate for times larger than about fdamp- 

corresponding LO equations are obtained by neglecting all contributions of order 
(d Xl ‘d PiJ ) in ©■ The LO equations are typically used to obtain kinetic equations 
with the further assumption of a quasi-particle picture. 

Fig. El shows the time evolution of the on-shell derivative 2 ujdxoF. The thick 
curve represents the “full” result, which is obtained from solving the evolution 
equations 0 and © for the three-loop 2PI effective action. For comparison, we 
evaluate the same quantity using the LO gradient expansion according to Eq. ©• 
For this we evaluate the RHS of the LO terms in Eq. © using the full result for F. 
If the gradient expansion to lowest order is correct, then both results have to agree. 
Indeed, one finds that the curves for the LO (dashed) and the full result indeed agree 
at sufficiently late times. Taking into account the NLO contributions the agreement 
can be improved substantially. The NLO curve (dotted) in Fig.EJapproaches the full 
result rather closely, however, they only agree after some characteristic time. The 
latter is determined by the time scale ~ fdamp for the effective loss of details about 
the initial conditions. This can be observed, e.g., from the decay of the unequal-time 
two-point function shown in the inset. The latter measures correlations at time t 
with the initial state and its decay on the time scale ~ fdamp coincides rather well 
with the time for the onset of validity of the NLO result. 

We have done a series of comparisons for various couplings in Ref. [ E| and 
conclude that for times sufficiently large compared to fdamp the gradient expansion 
seems to converge well. There are sizeable NLO corrections already for couplings of 
about A ~ 1/4. Nevertheless, one observes that the NLO result can get rather close 
to the full result for A < 1 [Ej. Times shorter than about fdamp seem clearly to be 
beyond the range of validity of transport equations even for weak couplings. This 
should make them unsuitable to discuss aspects of apparent early thermalization 
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and stresses the need to employ proper equations such as and © for initial-value 
problems. 


3. Beyond 2PI expansions: direct lattice simulations 


We have seen above that the techniques based on two-particle or higher irreducible 
generating functionals are crucial for our analytical understanding of nonequilibrium 
quantum fields. However, analytical approaches necessarily involve approximations 
such as a 2PI loop expansion. Nonequilibrium approximations are difficult to test 
for crucial questions of QCD, i.e. where strong interactions can play an important 
role. Here direct numerical simulations of the quantum held theory on a space- 
time lattice, i.e. without truncations, could boost our knowledge and trigger the 
development of further approximate analytical tools. 

Despite the importance of non-perturbative lattice simulation techniques in out- 
of-equilibrium quantum held theory, these are still in its infancies. This is in sharp 
contrast to well-established thermal equilibrium methods. Equilibrium calculations 
can typically be based on a Euclidean formulation, where the time variable is ana¬ 
lytically continued to imaginary values. By this the quantum theory is mapped onto 
a statistical mechanics problem, which can be simulated by importance sampling 
techniques. Nonequilibrium problems, however, are not amenable to a Euclidean 
formulation. Moreover, for real times standard importance sampling is not possible 
because of a non-positive definite probability measure. A very interesting recent de¬ 
velopment employs stochastic quantization techniques for real times, which do not 
require a probability distribution [ |Hj . Clearly, the possibility of direct simulations 
in nonequilibrium quantum held theory would mark a breakthrough not only for 
the description of QCD dynamics. 

In Ref. [ |U first lattice simulations of nonequilibrium quantum helds in Min¬ 
kowski space-time have been presented. For the example of a scalar held theory with 
quartic self-interaction this was used to compute the time evolution of correlation 
functions and characteristic time scales. Starting from a non-thermal initial state, 
the real-time quantum ensemble in 3+1 dimensions is constructed by a stochastic 
process in an additional (5th) “Langevin-time” using the reformulation of stochastic 
quantization for the Minkowskian path integral [El- In addition to the space-time 
variable x the held (f> depends on the Langevin-time parameter i) with 


d(j) 


■ SS[4>] M 

l ^ + V - 


( 8 ) 


where S denotes the classical action and 77 Gaussian or white noise. For a real 
quantum held theory the Langevin dynamics governs a complex 4> = (f)R + i(f>i, 
where the appearance of an imaginary part reflects the fact that in the quantum 
theory the held picks up a phase by evolving in time. 

The stochastic process © is associated to a real distribution P{cj) i 9) and 
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averages of observables A{<)>) are given as area integrals in the complex field plane: 

= f[d(j>R}[d(/)i}A((/) R + <j> 1 ','&) = /[d (j) R }A{(j) R )P eS ((t) R \•&) 

' /[d<fo][d0j]P(<fo,0/;??) _ /[d^fllPeff^i?;'!?) 

where the second equality defines the complex pseudo-distribution P e s{4>R', i?). The 
latter is indeed governed by the analytic continuation of the Fokker-Planck equation 
to real times, which admits the stationary solution [0 

lhnP eff (fe;^)~e iS N, (10) 

v —>oo 

Thus the approach can in principle be used for a Minkowskian theory, with “en¬ 
semble” averages calculated as averages along Langevin trajectories. 

Some properties seem to make the approach quite suitable for out-of-equilibrium 
calculations. Firstly, nonequilibrium requires specification of an initial state or 
density matrix. Therefore, the initial configuration is fixed which seems to stabilize 
the procedure. Moreover, the additional averaging over an initial density matrix can 
help to achieve fast convergence. Secondly, one typically has a good guess for the 
3 + 1 dimensional starting configurations of the Langevin updating procedure: In 
contrast to the quantum theory, the corresponding classical statistical field theory 
can be simulated using numerical integration and Monte Carlo techniques [ [5j. 
Using the nonequilibrium classical statistical solution as the starting configuration 
can improve convergence. Here the classical field configurations are obtained by 
numerically solving the classical field equation of motion and sampling over initial 
conditions, with nonzero field average and Gaussian fluctuations. It also provides 
a crucial check of the quantum result in some limiting cases: For sufficiently large 
macroscopic field or occupation numbers classical dynamics can provide a good 
approximation for the quantum evolution at not too late times [nn- 

Fig. 0 shows the time evolution for the connected part of the unequal-time 
correlator R e(f $(0, cc)<F(t, x)), which measures the correlation of the field at time 
t with the initial field. It gives important information about the characteristic time 
scale for the loss of details about the initial conditions, as explained in Sec. |T| 

One finds good convergence properties of the quantum simulations, which is a 
remarkable result. For given initial field configurations at time t = 0, very different 
starting configurations for the 3+1 dimensional space-time lattice converge to the 
same nonequilibrium dynamics for all t > 0. To obtain this one had to resolve the 
problem of possible unstable dynamics for the updating procedure [ |SJ . Though 
more or less formal proofs of equivalence of the stochastic approach and the path 
integral formulation have been given for Minkowski space-time, not much is known 
about the general convergence properties and its reliability. 

Two procedures can be employed for further tests of the algorithm. Firstly, one 
can compare to analytical approximations based on two-particle irreducible effective 
actions [El- Secondly, going to sufficiently late times one can compare to certain 
thermal equilibrium results from Euclidean simulations. 
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Fig. 5. The real part of (f $(0 ,x)<fr(t,x)) as a funtion of time t in units of the 
lattice spacing a. As starting configuration (i? = 0) the classical result is taken, and 
the Langevin updating (i9 > 0) incorporates quantum corrections [ jj£. 

The range of potential applications of first-principles simulations in nonequilib¬ 
rium quantum held theory is enormous. They may be used for out-of-equilibrium 
as well as Minkowskian equilibrium properties extracted at late times. Possible 
applications to QCD require implementation in a non-Abelian gauge theory, which 
is work in progress [EH 

We thank Ion-Olimpiu Stamatescu and Christof Wetterich for very fruitful col¬ 
laborations on recent related work. 
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